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Abstract 

We derive a method to study the phase diagram for high temperature super- 
conductors (HTCS). Our starting point is the Hubbard Hamiltonian with a 
weak attractive interaction to obtain the formation of bound pairs. We con- 
sider this attractive potential at different positions for different compounds 
accordingly to the experimental results of the coherence length. We then con- 
struct a wave function of the BCS type by a variational method using the 
Fourier transform of this extended Hubbard potential and then derive an en- 
ergy gap equation. This approach allows us to obtain the critical temperature 
as function of the doping concentration which gives very good agreement with 

the experimental phase diagrams of YBaCuO and La(Sr,Ba)CuO compounds. 
PACS: 74.20.Fg, 71.28+d, 74.70.tx 
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Most of the HTSC are quasi-two-dimensional insulators which become a metallic con- 
ductor and a superconductor below some critical temperature [|1],0. Another property very 
different from the usual superconductors is the very short coherence length of an electron 
(or hole) pair, ^ ^ 10 A. At present, there is no clear consensus about the origin of the 
mechanism of attraction but in several proposed models 0-0], the superconductor state is 
achieved from a hard core charged boson (formed by real space pairs) condensation in anal- 
ogy with the ^He problem. On the other hand, LDA calculations ^ indicates that the main 
features of the La2Cu04 band structure can be understood in terms of a two-dimensional 
tight-binding model. Another important point is that the properties of the normal metallic 
state are different than those of a common metal described by a free electron gas . As it is 
well known from the study of the Hubbard Hamiltonian , the on-site Coulomb correlations 
may explain the antiferromagnetism at low doping regime, the large magnetic fluctuations 
and the semiconductor-like properties of the metallic phase. On the other hand, in order 
to derive the formation of pairs and their binding energy, a weak attractive interaction Ui 
may be added to the Coulomb on-site repulsion Uq . The existence of bound states sug- 
gests that the normal ground state for many electrons in a tight-binding band may become 
unstable in the presence of these interactions Uq and Ui. 

In this letter, we shall use a BCS type wave function and a variational method to derive an 
energy gap or order parameter equation. We shall also use the intersite attractive potential 
Ui at different positions than the usual nearest-neighbor of the extended Hubbard model. 
These choices of intersite positions are directly determined from the experimental values of 
the coherence length and therefore, they depend on the speciflc compound to be studied. 
This procedure enable us to obtain the variation of Tc on the hole concentration which 
we compare with the experimental critical temperatures curves for the YBaCuO and the 
La(Sr,Ba)CuO compounds. Thus, let us start considering the two-dimensional extended 
Hubbard model on a square lattice 

(ij),a i (ij) 
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where t is the transfer integral (and the band width is 16t), Uq is the on-site Coulomb 
interaction and Ui is an intersite attraction and < ij > refers to nearest-neighbor pairs. 

Let us now study the two-electrons (or two-holes) problem. In this low-density limit, for 
s-wave pairs, an exact solution of the Schrodinger equation can be worked out in terms of 
the lattice Green's function. The biding energy A for a pair just below the bottom of the 
band is given by 

tr (Of I -Ui/2t-UoU,/m'-i 

tGooiSt + A) - _^u,/t-UoU,/2t^ + Uo/f 

To obtain A we make use of an expansion for Gqo in terms of an elliptic integral. In the 
limit of A <^ t, it was shown that |^ 
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tGoo{8t + A) ^ — [{1.38 + 0.2^A/t) - (0.25 + 0.125^A/t) ln(4A/t)j . (3) 

These equations allow us to study how A varies with Uo/t and Uq/Ui. In fact, Eq. (^ 
is valid only in the low density limit and it suggests the formation of electron (hole) pairs 
at the bottom of the band. A derivation of Eq. (|^), using real space methods 0, shows 
that these electron pairs have center of mass at rest and consequently, behave like Cooper 
pairs with momenta k and —k. After these preliminary considerations, let us focus on the 
many-body problem. For this purpose, we construct a trial wave function of the BCS type 
describing pair of electrons (or holes) 

k 

where |$o) is the empty band state and m| -|- f | = 1. 

Following the variational approach [jlO[, we must minimize the expression 



-/i($|iV|$), (5) 

which yields 

{m-m = 2i: + e ^^z^.^wn (6) 
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with = —4:t{cos{kxa) +cos{kya)) — fi and with V^^-being the interaction part of the potential 
that describes the transition of a pair from the state {k, ~k) to (/, — /). The minimization 
procedure follows exactly as the BCS theory for the free electron gas [jlO[ and we obtain the 
same type of T = energy gap equation 
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As we already mentioned, V^f is the Fourier transform of the potential of Eq. (|I|), which 
is approximately given by 

= [Vo^^ - 2^'^' cos(A;,. - Qa] [V^^'' - 2Vl^^ cos(A;, - /,)a] , (8) 
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Comparing with the Hamiltonian Eq. we identify Ui = \/VoVi and Uq = Vq and we 
use that Vq > Vi. Likewise BCS theory, we assume the gap to have the same functional 
form of the potential, namely, = A(0)(cos(/c^a) + cos{kya))/2. This functional form was 



also previously deduced for the RVB particle-particle energy gap |TT[| using self-consistent 



field methods ||T0[. Although this potential-energy gap relation is the same used in the BCS 



theory, it is worth to stress that the above potential is very different than the constant and 
isotropic mean potential of the original BCS theory and it has the correlations built on. We 
take /i equal to the hole maximum energy and we will assume that it grows linearly with the 
concentration of holes. Although this dependence is an ansatz here, it has been previously 
derived by auxiliary-bosons mean field theory [|ll|,|l3 . 



We calculate the probability of finding a hole pair, that is, the condensation amplitude 
Fj: = u^vj:. It has a maximum at fc^ and drops very rapidly with \k\ > kM, that is, the 
pair formation instability occurs mostly at the Fermi surface. This is the same result of 
the free-electron BCS theory, despite the fact that our potential (Eq. (^)) acts on all the 
first Brillouin zone. This is mainly due to the small values of the Fermi surface and the 



very short coherence length which indicates that all carriers can be involved in pairing [Q. 
According to these considerations and taking A; = in Eq. (|^), we obtain 

1 = / / ^ ^ -j^ dadp, (10) 

where f{a,P) = cos a + cos/3 and a = kxa and P = kya. and Pm are the maximum 
T = occupied values (like a Fermi momentum) that depend on the density of holes or, 
as more currently used, on the number of holes per Cu atoms, x. Then we can show that 
aM = Pm ~ arcos{l—x). The integrations are performed up to a a and P^, which are chosen 
at values where the condensation amplitude becomes very small, namely, ~ 0.01. This is 
usually attained for > 6A(0) and larger values of k do not modify the final results. The 
two-dimensional integrals are done by an elementary Simpson's rule algorithm [|l^. Thus, 
for given values of Vi/Vq and Vo/t, there are corresponding values of A(0) and clearly the 
relation among them depends on x. Furthermore, there may be only real solutions at a 
certain range of x. 

We have derived an expression for the T = energy gap (Eq. (PH])). For T 7^ 0, the 
excitations with their respective probability must be taken into account. The derivation 
of a self-consistently temperature-dependent gap equation is analogous to that which leads 
to Eq. (|10|). At this point, we again follow the BCS approach |jlO| and assume that A(r) 



vanishes at the critical temperature Tc, which yields the following equation 

1 ra^ rP^ (Fo - 4 {ViV.f^ f{a, P)) f{a, P) tanh 



2vr2 7o Jo At{f{a,P)-f{aM,PM)) ^ ^' ^ ^ 

where we again integrate up to oa and Pa which except for x close to one, is almost entirely 
within the lower Hubbard band. It is a well known result that the correlation Vq can split 
the half-filled conduction band into two, with the upper band being empty and the lower 



band being filled |T^ . Accordingly the properties of our equations are also dominated by the 
lower band. We should also point out that in the derivation of Eq.(ll), we used the entropy 
form of a non- interacting fermion system which is an approximation which is justified in the 
low energy gap limit, that is, for A(0) ^ 16t. 



Now we compare Eq. ([T0|) and Eq. (|TlD in order to relate Vi/Vq, Vo/t, A(0) and Tc. The 
first important thing to notice is that the different coherence lengths for different compounds 
suggests that the intersite attractive pair potential should be placed at a different site than 
the original nearest neighbor site of the extended Hubbard Hamiltonian. In other words, the 
minimum of the real attractive potential depends on the type of HTCS. This can be easily 
included in our equations by a change in the lattice parameter only in the potential (Eq. 
(I) and expressions. Thus, for the lanthanum compounds, since the lattice parameter 
is approximately AA and ^ ~ 35A P] , we expect that the minimum of the attractive pair 
potential occurs at a distance of the order of 8 — 9 lattice parameters. Similarly, with a 
coherent length ^ ^ 15 A p|, we expect the attractive pair potential minimum to be at 3 — 4 
lattice parameters for the yttrium compounds. Thus, in order to obtain values pertinent to 
the experimental results for the La2-xSrj.Cu04 compounds, we use the attractive potential 
at the 6*'*-neighbor position and with Vq/Vi = 8^. We find a real solution only for values 
of X G [0.04,0.35] and with a maximum value for Tc at x = 0.16. Thus, the position of 
the maximum is fixed and the absolute values of Tc depend only on A(0). Therefore, to 
compare with the measurements, we use 16t = 2ev and we choose A(0) = 80i^' in order 
to obtain Tc = 35K at x = 0.16. In Fig.l, we plotted our values for Tc as function of x and 
we see that they provide an excellent fit for the experimental data points which were taken 
from Refs. Furthermore, we obtain values for greater than 2.3 (around x = 0.16) 



which agrees with earlier measurements [0 which yielded values larger than 2.5. These 
results are considerably larger than the BCS value of 1.75. 

For YBa2Cu30Q+x we place the attractive potential at the next neighbor position. With 
Vo/^i = 5^, we find that there are solutions only for x > 0.4 and the maximum Tc occur 
at X = 0.96, which is close to the measured maximum Tc at x = 0.93. Now we choose 
A(0) = 1807^^ which gives Tc = 95K at x = 0.96. Our results are plotted in Fig.2 with 
also some experimental data points taken from Ref. We notice that the agreement in 
this case is only qualitative and fail to reproduce the change in concavity around x = 0.7. 
A possibility is that this structure is due to inter-plane effects in the compounds richer in 



oxygen and, therefore, it should not be reproduced by our two-dimensional treatment. A 
difficulty with our calculations is the variation of the values of Vo/t (and Vi/t) with x. For 
the lanthanum compounds, as a result of Eqs.(lO) and (11), we obtain Vo/t = 11 for x = 0.16 
and at the onset of superconductivity, we find Vo/t ^ 11 (at x = 0.05). Although they are in 
the same order of magnitude of the expected values of Vo/t = 16 such variation of the 
coupling with the concentration is artificial. It is probably due to the linear approximation 
used for the chemical potential. 

In summary, we demonstrated that a variational procedure with a BCS type wave func- 
tion in connection with a two-dimensional Hubbard Hamiltonian and with an attractive 
interaction is suitable to determine the superconducting state for layered HTSC This ap- 
proach allows us to derive the energy gap equation in fc-space, which is used to relate 
with X. The results yield an excellent ffiting for the curves of the critical temperature as 
function of doping when we take an effective position for the attractive potential at a dis- 
tance consistent with the measured values for the coherence length. The different positions 
for the attractive potential are probably due to the coupling to a bosonic field (phonons, 
excitons, plasmons, etc....). Such field is polarized along the electronic motion and induces 
an effective short range attraction (renormalized from its bare value) which varies for differ- 
ent compounds. Other experimental values, like the ratio A/Tc, are well reproduced by our 
method. Thus, we conclude that the method derived here can be successfully applied to the 
two most studied compounds, even though, they have transitions at very different ranges of 
T and x. 
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Mauro Doria. He also thanks CNPq and FINEP for partial financial support. 
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FIGURES 

FIG. 1. The critical temperature for La2-xSrxCu04 as function of the hole concentration x. 

The solid line is plotted from our calculations and the black dots are experimental results taken 
from Refs.4 and 5 . 

FIG. 2. The critical temperature for YBa2CusOQ+x as function of the hole concentration x. 
The solid line is plotted from our calculations and the black dots are experimental results taken 
from a plot in Ref.2 . 
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